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^ ; ABSTRACT 

, We perform a joint determination of the distance-redshift relation and cosmic expan- 

' sion rate at redshifts z = 0.44, 0.6 and 0.73 by combining measurements of the baryon 

I acoustic peak and Alcock-Paczynski distortion from galaxy clustering in the WiggleZ 

, Dark Energy Survey, using a large ensemble of mock catalogues to calculate the covari- 

• ance between the measurements. We find that Da (-2) = (1205 ± 114, 1380 ± 95, 1534 ± 

• '-j ; 107) Mpc and H{z) = (82.6 ± 7.8, 87.9 ±6.1, 97.3 ± 7.0) km s^^ Mpc^i at these three 

p% . redshifts. Further combining our results with other baryon acoustic oscillation and dis- 

tant supernovae datasets, we use a Monte Carlo Markov Chain technique to determine 
the evolution of the Hubble parameter H{z) as a stepwise function in 9 redshift bins 
of width Az = 0.1, also marginalizing over the spatial curvature. Our measurements 
of H{z), which have precision better than 7% in most redshift bins, are consistent 
with the expansion history predicted by a cosmological-constant dark-energy model, 
in which the expansion rate accelerates at redshift z < 0.7. 
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1 INTRODUCTION 

One of the fundamental goals of observational cosmology is 
to determine the expansion rate of the Universe as a function 
of redshift. Measurements of the expansion history, which 
can be described by the evolution of the Hubble parameter 
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H{z) = (1 + 2) da/dt with redshift z, where a{t) is the cos- 
mic scale factor at time t, provide one of the most important 
observational tests of the cosmological models which char- 
acterize the different components of the Universe and their 
evolution with time. In particular, a paramount problem in 
cosmology is to understand the physical significance of the 
"dark energy" which appears to dominate the cosmic energy 
density today, as described by the phenomenology of the 
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standard cosmological-constant cold-dark-matter (ACDM) 
model. 

A number of powerful tools to measure the cosmic ex- 
pansion history beyond the local Universe have been de- 
veloped in recent decades. Foremost amongst these probes 
is the use of distant Type la supernovae (SNe) as standard 
candles (e.g. Riess et al. 1998, Perlmutter et al. 1999, Aman- 
ullah ot al. 2010, Conley ct al. 2011, Suzuki ct al. 2011). 
The apparent peak magnitude of these supernovae, follow- 
ing certain corrections based on the light-curve shape which 
decrease the observed scatter in the peak brightness, yield 
a relative luminosity distance as a function of redshift, i.e. 
Dl{z) Ho/c where Dl{z) is the luminosity distance, Hq is 
the local Hubble parameter and c is the speed of light. 

Although such measurements accurately trace the 
shape of the distanco-rodshift "Hubble diagram" in the red- 
shift range z < 1, a, number of qualifications must be men- 
tioned. Firstly, the expansion history H{z) is not directly 
measured by supernovae but must be determined as a deriva- 
tive of the noisy luminosity distances (Wang & Tegmark 
2005, SoUerman et al. 2009, Shaficloo & Clarkson 2010). 
Secondly, obtaining the expansion rate from the luminos- 
ity distance requires an additional assumption about spatial 
curvature, which influences the geodesies followed by pho- 
tons. Thirdly, despite the impressive and thorough treat- 
ment in recent supernovae analyses of the systematic errors 
which could bias cosmological fits, these systematics now 
limit the utility of these datasets. 

Large galaxy surveys offer a complementary route for 
mapping cosmic distances and expansion, using two prin- 
cipal techniques. Firstly, the largo-scale clustering pattern 
of galajcies contains the signature of baryon acoustic oscil- 
lations (BAOs), a preferred length scale imprinted in the 
distribution of photons and baryons by the propagation of 
sound waves in the relativistic plasma of the early Universe. 
This length scale, the sound horizon at the baryon drag 
epoch rs{zd), may be accurately calibrated by observations 
of the Cosmic Microwave Background (CMB) radiation and 
applied as a cosmological standard ruler (Blake & Glaze- 
brook 2003, Seo & Eisenstein 2003). Some applications to 
galaxy datasets are presented by Eisenstein et al. (2005), 
Percival et al. (2010), Beutler et al. (2011), Blake et al. 
(2011b), Padmanabhan et al. (2012) and Anderson et al. 
(2012). 

Given a sufficiently-large galaxy survey at a redshift z, 
the preferred scale may be detected in both the tangential 
direction on the sky as an enhancement in the number of 
galaxy pairs with a given angular separation A^, and in the 
radial direction as an excess of pairs with redshift separation 
Az. If these two signals can be simultaneously extracted by 
measuring galaxy clustering in tangential and radial bins, 
they respectively carry information about the angular di- 
ameter distance Da{z) — Dl{z)/{1 + z)'^ and the Hubble 
expansion rate at the survey redshift in units of the standard 
ruler, r,{zd)/[{l + z)Da{z)] ~ AO and r,{zd)H{z) ~ cAz. 
If the galaxy survey only permits the baryon acoustic peak 
to be detected in the angle-averaged galaxy clustering pat- 
tern, then an effective "dilation scale" distance is measured 
which consists of two parts Da(2) and one part 1/H{z): 
Dv{z) = [{l+zfDA{zfcz/H{z)]^^^ (Eisenstein et al. 2005, 
Padmanabhan & White 2008). 

The second technique through which large-scale struc- 



ture surveys permit measurement of geometrical distances 
is the "Alcock-Paczynski (AP) test" (Alcock & Paczynski 
1979). The AP test probes the cosmological model by com- 
paring the observed tangential and radial dimensions of ob- 
jects which are assumed to be isotropic in the correct choice 
of model. It can be applied to the 2-point statistics of galaxy 
clustering if redshift-space distortions, the principal addi- 
tional source of anisotropy, can be successfully modelled 
(Ballinger, Peacock & Heavens 1996, Matsubara & Suto 
1996, Simpson & Peacock 2010). By equating radial and 
tangential physical scales, independently of any underlying 
standard ruler, the observable Az/ Ad ~ (1+z)Da{z)H{z)/c 
may be determined (Outram et al. 2004, Marinoni & Buzzi 
2010, Blake et al. 2011c). 

Therefore, using a combination of BAO and/or AP 
measurements, large-scale galaxy surveys can supply in- 
dependent measurements of the distancc-rcdshift relation 
Da{z) and expansion history H{z). We note that the Hub- 
ble expansion rate as a function of redshift may also be in- 
ferred from the relative ages of passively-evolving galaxies 
(Jimenez & Loeb 2002, Stern et al. 2010, Carson & Nichol 
2010, Moresco et al. 2012). This is a promising technique 
albeit subject to assumptions about the stellar populations 
of these galaxies, in particular about galaxy metallicity at 
high redshift. 

We focus here on distance measurements using the Wig- 
gleZ Dark Energy Survey (Drinkwater et al. 2010), which 
was designed to extend the study of large-scale structure 
over large cosmic volumes to redshifts z > 0.5. The study 
presented here builds upon two existing distance-scale mea- 
surements using the WigglcZ dataset: Blake et al. (2011b) 
reported the measurement of the angle- averaged baryon 
acoustic peak at redshifts z = (0.44, 0.6, 0.73), and Blake et 
al. (2011c) applied the Alcock-Paczynski test to the 2D clus- 
tering power spectrum. In this latter study we combined the 
AP fits with SNe data to estimate the Hubble parameter rel- 
ative to its local value, H{z)/Ho. We here extend these anal- 
yses by combining the WiggleZ measurements of D\/H and 
DaH , including a calculation of the covariance of the statis- 
tics, to extract measurements of Da{z) and H{z) as a func- 
tion of redshift, in absolute units independent of the value 
of Hq, based solely on WiggleZ Survey data (and a sound- 
horizon calibration). F\irthermore, by combining these mea- 
surements with SNe and other galaxy datasets, we constrain 
the cosmic expansion rate H{z) as a stepwise function in the 
redshift range z < 0.9, and fit a variety of cosmological mod- 
els to the results. 



2 DATA 

2.1 WiggleZ Survey clustering measurements 

The WiggleZ Dark Energy Survey (Drinkwater ct al. 2010) 
is a large-scale redshift survey of bright emission-line galax- 
ies which was carried out at the Anglo- Australian Telescope 
between August 2006 and January 2011. The galaxy sam- 
ple utilized by this study is drawn from the final set of ob- 
servations covering about 800 deg^ of sky in six regions, 
including a total oi N = 158,741 galaxies in the redshift 
range 0.2 < z < 1.0, and is the same dataset as used for 
the analysis of the baryon acoustic peak by Blake et al. 
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Figure 1. The joint likelihood of the Alcock-Paczynski scale distortion parameter F(z) = (1+z)Da{z)H{z)/c and the normalized growth 
rate quantified by f cr^{z), obtained from fits to the 2D galaxy power spectra of the WiggleZ Dark Energy Survey in three overlapping 
redshift slices 0.2 < z < 0.6, 0.4 < z < 0.8 and 0.6 < z < 1.0. This Figure was produced by marginalizing over the linear bias factor 
. The probability density is plotted as contours enclosing 68.27% and 95.45% of the total likelihood. The solid circles indicate the 
parameter values in a fiducial flat ACDM cosmological model with parameters f2ni = 0.27, trg = 0.8. 



(2011b). Our study is based on measurements of the angle- 
averaged galaxy correlation function and 2D galaxy power 
spectrum in tangential and radial Fourier bins in three over- 
lapping redshift ranges 0.2 < z < 0.6, 0.4 < z < 0.8 and 
0.6 < z < 1.0. The effective redshifts of the measurements 
in these three redshift slices are z^a = (0.44, 0.6, 0.73) (Blake 
et al. 2011b). We use overlapping, wide redshift ranges in or- 
der to ensure a detection of the baryon acoustic peak in each 
redshift shce (following Percival et al. 2010) and to provide 
the best mapping of the distance-redshift relation. We ac- 
count for the correlations between the measurements when 
fitting models. 

2.2 Fitting the baryon acoustic peak 

Blake et al. (2011b) presented our analysis of the WiggleZ 
galaxy correlation function to map the baryon acoustic peak 
in these three redshift ranges. The correlation functions con- 
tain evidence for the baryon acoustic peak in each redshift 
slice. Our model for fitting the correlation function measure- 
ments to determine the standard-ruler distance is described 
in Section 3.1 of Blake et al. (2011b). Our results are most 
cleanly expressed as a measurement of the acoustic param- 
eter, 

_ 100^)v(^)^/^O? 
A[z) - - , (1) 

at the effective redshift of the sample. We marginalized over 
the shape of the clustering pattern (parameterized by the 
physical matter density f2mft^), non-linear damping of the 
acoustic peak, and galaxy bias. The results for A[z) are 
listed in Table 1, reproduced from Blake et al. (2011b). 

2.3 Fitting the 2D power spectrum 

Blake et al. (2011c) presented an analysis of the 2D WiggleZ 
galaxy power spectrum, where modes are binned by Fourier 
wavenumber and angle to the line-of-sight. We repeated 
these measurements for the new choice of redshift bins con- 
sistent with the correlation function analysis. Our model 
for fitting the 2D power spectrum to extract the Alcock- 
Paczynski distortion is described in Section 3.2 of Blake et 
al. (2011c). At each redshift we obtain a measurement of the 



distortion parameter F(z) = (1 -f z)Da{z)H{z)/c and the 
normalized growth rate / as, which quantifies the amplitude 
of redshift-space distortions in terms of the growth rate / 
and amplitude of matter fluctuations ag . We also marginal- 
ized over a linear bias factor which we do not require to 
be identical to the bias in the correlation function model. 
This is a conservative choice which reflects the possibility 
that the amplitude of these two statistics due to galaxy bias 
and redshift-space distortions may be scale-dependent. Fig- 
ure 1 displays the joint likelihoods of F and / erg fitted in 
each of the new redshift slices. We note that although there 
is a strong correlation between the parameters, both may 
be successfully determined, and the results are collected in 
Table 1. 



3 JOINT FITS FOR THE EXPANSION AND 
GROWTH HISTORY 

3.1 Covariances between fitted parameters 

Given that our measurements of the baryon acoustic scale 
and tangential/radial clustering anisotropy have taken place 
in overlapping redshift slices, using clustering statistics 
which are potentially correlated by common cosmic vari- 
ance, it is important to determine the covariances between 
the measurements of {A, F, fas) within and between redshift 
slices. This is achieved by repeating the different parameter 
fits for each of a series of lognormal realizations, and using 
the statistical ensemble to determine the various correla- 
tion coefficients. We generated 400 lognormal realizations 
for each WiggleZ survey region and redshift slice (i.e., 2400 
realizations for each redshift slice, or 7200 in total) using 
the methods described by Blake et al. (2011a). Lognormal 
realizations provide a reasonably accurate galaxy clustering 
model for the linear and quasi-linear scales which are im- 
portant for modelling the large-scale clustering pattern. 

Figure 2 displays the correlations between single param- 
eters fit to different pairs of redshift slices. In each panel, 
the small dots represent the best-fitting pairs of parameter 
values in the redshift slices for the 400 lognormal realiza- 
tions. The red ellipses are 2D Gaussians representing the 
covariance between the fitted parameters, and the solid red 
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Table 1 . Results of cosmological model fits to the galaxy correlation functions and power spectra measured in three overlapping redshift 
slices of the WiggleZ Dark Energy Survey. The acoustic scale parameter A{z) is obtained from the fit to the baryon acoustic peak in 
the correlation function (marginalizing over the matter density QjnhP, the damping parameter <t„ and the galaxy bias factor) and the 
parameters F(z) and f <ts(z) are measured using the 2D power spectra (marginalizing over a galaxy bias which is not assumed to be 
identical to its value in the correlation function fit). The angular diameter distance Dj^{z) and Hubble expansion rate H{z) are derived 
from the measurements of A.(z) and F(z) assuming a CMB-motivated prior in Qm/i^ in order to calibrate the standard ruler. 



Redshift slice 0.2 < z < 0.6 0.4 < z < 0.8 0.6 < z < 1.0 

Effective redshift z 0.44 0.60 0.73 

A(z) = lOODv (z) VTwP'/cz 0.474 ± 0.034 0.442 ± 0.020 0.424 ± 0.021 

F{z) s (1 + z)Da{z)H{z)/c 0.482 ± 0.049 0.650 ± 0.053 0.865 ± 0.073 

/o-g(z) 0.413 ±0.080 0.390 ±0.063 0.437 ± 0.072 

Da{z) [Mpc] 1204.9 ±113.6 1380.1 ±94.8 1533.7 ±106.8 

i/(z) [km s-i Mpc-i] 82.6 ±7.8 87.9 ±6.1 97.3 ± 7.0 



circle is the input fiducial model used to generate the log- 
normal realizations. The correlation coefficient r is quoted 
in the bottom left-hand corner of each panel, and is consis- 
tent with zero in the second column when non-overlapping 
redshift slices are used. 

Figure 3 shows the correlations between pairs of diflier- 
ent parameters fit in the same redshift slice, using the same 
presentation format as Figure 2. The strongest covariance 
is measured between the Alcock-Paczynski distortion F and 
growth rate / as from redshift-space distortions, with corre- 
lation coefficients r ~ 0.8. The measurements of the baryon 
acoustic peak "monopole" parameter A are correlated with 
each of the "quadrupole" parameters {F, fas) at a lower 
level r ~ 0.2. The full 9x9 covariance matrix for the pa- 
rameters is listed in Table 2. 

Figure 4 plots the ID distributions of best-fitting pa- 
rameters for the second redshift slice, demonstrating that 
this is well-described by the multivariate Gaussian model 
and does not contain significant wings that might cause the 
confidence regions to be under-estimated in subsequent cos- 
mological parameter fits. The distributions for the other red- 
shift slices are similar. 



(0.44,0.60,0.73) are Da{z) = (1205±114, 1380±95, 1534± 
107) Mpc and H{z) = (82.6 ± 7.8,87.9 ± 6.1,97.3 ± 7.0) 
km s~^ Mpc~^. A steady increase in the value of H{z) 
with z is consistent with accelerating expansion given that 
dH/dt = —H^ [1 + (a/a)] is negative when a > 0. These 
measurements of Da{z) and H{z) are listed in Table 1 
along with the marginalized measurements of the normal- 
ized growth rate. The fractional accuracies with which the 
parameters are measured arc 7 — 9% for Da and H, and 
16 — 20% for f as- We note that readers wishing to include 
our dataset in cosmological parameter fits should use the 
raw measurements of {A,F,fas) given in Table 1, together 
with the covariance matrix listed in Table 2, rather than 
these derived values of Da and H. 



4 COSMOLOGICAL MODEL FITS 

We now use our joint WiggleZ measurements of the baryon 
acoustic peak and Alcock-Paczynski distortions to place 
constraints on parametric and non-parametric cosmological 
models, both alone and in combination with other datasets. 



3.2 Determination of Da{z) and H{z) 

Using the joint measurements of the Alcock-Paczynski dis- 
tortion parameter F oc DaH and acoustic parameter A oc 
{D\/HY'^ in each redshift shce, we can break the degener- 
acy between the angular-diameter distance Da{z) and Hub- 
ble parameter H(z). We fit for Da and H in each redshift 
slice using these measurements and their covariance. We also 
marginalize over the physical matter density fimft^, which 
appears in Equation 1 for A{z), using a Gaussian prior with 
mean 0.1345 and width 0.0055. This prior is motivated by 
fits to the CMB (Komatsu ct al. 2009) and is independent 
of the low-redshift expansion history under certain general 
assumptions, which are listed in Section 5.4.1 of Komatsu 
et al. (2009). 

The joint likelihood of Da{z) and H{z) in each of the 
three redshift slices is displayed in Figure 5, where the solid 
line represents the joint variation of these parameters with 
redshift in a fiducial cosmological model fim = 0.27 and 
h = 0.71, and the solid circles superimposed on the line in- 
dicate the model prediction for the three analyzed redshift 
slices. The marginalized values of Da and H at redshifts z = 



4.1 Other cosmological datasets 

Together with the new WiggleZ results described in this 
study, we add BAO distance measurements obtained from 
the 6-degree Field Galaxy Survey (6dFGS; Beutler et al. 
2011) and by applying "reconstruction" to the sample of 
Sloan Digital Sky Survey (SDSS) Luminous Red Galaxies 
(Padmanabhan ct al. 2012). We also include the joint BAO 
and AP measurements recently reported by the Baryon Os- 
cillation Spectroscopic Survey (BOSS, Reid et al. 2012). 

We additionally use the "Union 2" compilation of su- 
pernovae data by AmanuUah ct al. (2010), as obtained from 
the website http://supernova.lbl.gov/Union. This com- 
pilation of 557 supernovac includes data from Hanmy ct al. 
(1996), Riess et al. (1999, 2007), Astier et al. (2006), Jha ct 
al. (2006), Wood-Vasey et al. (2007), Holtzman et al. (2008), 
Hickcn ct al. (2009) and Kcsslcr ct al. (2009). When fitting 
cosmological models to this SNc dataset we used the full 
covariance matrix of these measurements including system- 
atic errors, as reported by AmanuUah et al. (2010). We also 
performed an analytic marginalization over the unknown ab- 
solute normalization (Goliath et al. 2001, Bridle et al. 2002). 
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Figure 2. Correlations of the parameters A{z), F(z) and f ag{z) when each parameter is fitted to pairs of the three overlapping WiggleZ 
rcdshift slices. The Alcock-Paczynski distortion parameter F is plotted relative to its value in the fiducial cosmology, Fa^. Each small dot 
represents the best-fitting values of the parameters using the correlation functions and power spectra measured from 400 independent 
lognormal realizations. The red ellipses represent the derived covariances between the parameter fits, and the solid red circle is the input 
fiducial model of the lognormal realizations. The correlation coefficient r is quoted in the bottom left-hand corner of each panel, and is 
consistent with zero in the second column when non-overlapping redshift slices are used. 



Table 2. This Table lists the values of 10^£Z, where d is the covariance matrix of measurements from the WiggleZ Survey data of 

the acoustic parameter A{z), the Alcock-Paczynski distortion parameter F{z) and normalized growth rate fas{z), where each pa- 
rameter is measured in three overlapping redshift slices (2:1,22,^3) with effective redshifts z^ff = 0.44, 0.6 and 0.73, respectively, where 
zi = [0.2,0.6], 22 = [0.4,0.8] and 23 = [0.6,1.0]. The data vector is ordered such that y^j^^ = {Ai, A2, A3, Fi, F2, F3, faa,ijas,2, ft^a.a) = 
(0.474,0.442,0.424,0.482,0.650,0.865,0.413,0.390,0.437). The chi-squared statistic for any cosmological model vector Y^^^ can be ob- 
tained via the matrix multiplication = Qlobs ~ Xmod)^C^~^ (^obs ~ Xmod)- ^^'^ matrix is symmetric; we just quote the upper 
diagonal. 



Parameter 


A(2i) 


A{Z2) 


A{Z3) 


^(^i) 


Fiz2) 


F{Z3) 




/o'8(22) 




A(zi) 


1.156 


0.211 


0.000 


0.400 


0.234 


0.000 


0.598 


0.129 


0.000 


A{z2) 




0.400 


0.189 


0.118 


0.276 


0.336 


0.080 


0.227 


0.230 


A{zz) 






0.441 


0.000 


0.167 


0.399 


0.000 


0.146 


0.333 


n-^i) 








2.401 


1.350 


0.000 


2.862 


1.080 


0.000 


F(Z2) 










2.809 


1.934 


1.611 


2.471 


1.641 


F{Z3) 












5.329 


0.000 


1.978 


4.468 
















6.400 


2.570 


0.000 


fc8{z2) 
















3.969 


2.540 


/ct8(z3) 


















5.184 



© 0000 RAS, MNRAS 000, 000-000 



6 Blake et al. 




0.4 0.5 
A(0.2<z<0.6) 




0.4 0.5 
A(0.2<z<0.6) 



^ 6 

d m 
V ■ 
N o 
V 

d 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — TTT- 




.r.= p/3 , , 

-I I I I I I I I I I I I I I I I I I l_ 



0.8 0.9 1 1.1 1.2 
F/Fpj(0.2<z<0.6) 



CO ^ 

d 

V 

N 
V 

■>i- ^ 
d 

V ' 

■D 



A vs F/F,i, 



r = 0.26 , 

_l I I I I I I I I ■_ 



^.3 0.4 0.5 

A(0.4<z<0.8) 



^ d 
oq 

d in 
V • 
N o 
V 

GO 

b 

d 



— I 1 1 1 1 1 1 1 1 1 1 r- 

A VS fa» 



r.=. 0.18 

I I I I I I I I I I I L. 



0.3 0.4 0.5 

A(0.4<z<0.8) 



j-i — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

oo 




V ' - 
N o F 

V 



o . r = 0,74 

Lj I I I I I I I I I I I I I I I I I ■_ 



0.8 0.9 1 1.1 1.2 
F/Ffl/0.4<z<0.8) 




0.3 0.35 0.4 0.45 0.5 
A(0.6<z<1.0) 



J I I I I I I I I I I I I I I I I I I I I I I I 



^ d 
q 

N o 
V 

S d 

CO 

b 

d 



(O 
^ d 



N o 
V 



A vs f(Tg . 




■ r = 0.22 ■ 

' I ■ ■ ■ ■ I I I I ■ ■ I I I I I ■ I I I I 



0.3 0.35 0.4 0.45 0.5 
A(0.6<z<1.0) 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

F/Fflj vs ftTg : . 




r = 0,85 

-I I I I I I I I I I I I I I I I I I ■_ 



0.8 0.9 1 1.1 1.2 
F/Fflj(0.6<z<1.0) 



Figure 3. Correlations between different pairs of the parameters A{z), F{z) and f a^(z) fitted to tfie tfiree WiggleZ redshift slices. 
The Alcock-Paczynski distortion parameter F is plotted relative to its value in the fiducial cosmology, -Ffid- Each small dot represents 
the best-fitting values of the parameters using the correlation functions and power spectra measured from 400 independent lognormal 
realizations. The red ellipses represent the derived covariances between the measurements, and the solid red circle is the input fiducial 
model of the lognormal realizations. The correlation coefficient r is quoted in the bottom left-hand corner of each panel. Although 
the strongest correlation is obtained between F{z) and f C7^(z), weaker but non-zero correlations are measured between both of these 
parameters and A[z). 
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Figure 4. The ID distribution of the best-fitting parameters for 400 lognormal realizations for the second redshift slice, compared to 
the adopted multivariate Gaussian model for the covariance. The full distribution does not contain significant wings that might cause 
confidence regions to be under-estimated in subsequent cosmological parameter fits. 
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Finally, in some fits we include CMB data using the 
Wilkinson Microwave Anisotropy Probe (WMAP) "distance 
priors" (Komatsu et al. 2009) using the 7-year WMAP re- 
sults (Komatsu et al. 2011). The distance priors quantify 
the complete CMB likelihood via a 3-parameter covariance 
matrix for the acoustic index £a, the shift parameter TZ and 
the redshift of recombination Zt, as given in Table 10 of 
Komatsu et al. (2011). When deriving these quantities we 
assumed a physical baryon density flhh^ = 0.0226, a CMB 
temperature Tcmb = 2.7257^ and a number of relativistic 
degrees of freedom A'^off ~ 3.04. 



4.2 loCDM fits to WiggleZ measurements 

As an initial analysis we fitted a flat wCDM cosmological 
model to these datasets in which spatial curvature is fixed 
at Jlk = but the equation-of-state w of dark energy is var- 
ied as a free parameter. We fitted for the three parameters 
(Qm, flmh^ , w) using flat, wide priors which extend well be- 
yond the regions of high likelihood and have no effect on the 
cosmological fits. We only use the joint WiggleZ measure- 
ments of A{z) and F{z) in these fits, not the growth rate 
data. The extra complexity in the normalization of the clus- 
tering pattern required to fit f as{z) is analyzed by Parkin- 
son et al. (in prep.). 

Figure 6 displays the joint likelihood of the parame- 
ters {Q,m, w) marginalizing over ilmft^, comparing the effects 
of adding different datasets to the WMAP distance priors. 
The combination with the joint WiggleZ measurements of 
A{z) and -^(2) is illustrated by the (black) solid contours, 
with the (red) dashed contours showing the improvement 
compared to only using the WiggleZ measurements of A[z). 
The results are consistent with, albeit with a significantly 
lower accuracy than, parameter measurements based on the 
combination of the WMAP distance priors with all Wig- 
gleZ-f-BOSS AP-fBAO data, and all SNe data, which are 
represented by the (blue) dash-dotted contours and (ma- 
genta) dotted contours, respectively. 
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Figure 5. The joint likelihood of fits of Da{z) and H{z) to the 
baryon acoustic peak and Alcock-Paczynski distortions in each of 
three overlapping WiggleZ redshift slices. The two contour levels 
in each case enclose regions containing 68.27% and 95.45% of the 
total likelihood. A flat ACDM model prediction for cosmological 
parameters Sim = 0.27 and h = 0.71 is plotted as the solid line, 
with circles representing the effective redshifts of the three data 
slices. 
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4.3 ll{z^ fits in bins 

We performed a model-independent determination of the 
cosmic expansion history by carrying out a Monte Carlo 
Markov Chain (MCMC) fit for the Hubble parameter H{z) 
as a stepwise function in narrow redshift bins, where we 
identify the value in the first bin as H(zi) — Ho ~ 100 h 
km s~^ Mpc~^. We also included the spatial curvature 
as a free parameter in our fit, such that we derived the an- 
gular diameter distance Da(z) from the co- moving radial 
co-ordinate r{z) as: 



DA{z)^r{z)/{l 

Da{z) = i?curv Sinh (r(2:)/7?curv)/(l 

Da{z) = RcuTv sin(r(2:)/7?curv)/(l 



J7k = 
J7k > 
J7k < (2) 



c/ (Ho\/|f2k|). If the range of the zth redshift 
and the redshift at which we 



< Z < Zi, 



where Rcurv 

bin is 2i,min 

are evaluating a distance lies in the nth bin, the co-moving 
radial co-ordinate r{z) is deduced from the stepwise H{z) 



Figure 6. The joint probability for parameters Qm and w fit- 
ted to different datasets added to the WMAP distance priors, 
marginalizing over Q^^h^ and assuming f2k = 0. The CMB data 
is combined in turn with the WiggleZ acoustic scale and Alcock- 
Paczynski distortion measurements (with appropriate covariance; 
black solid contours), the WiggleZ acoustic scale measurements 
alone (red dashed contours), all BAO measurements (blue dash- 
dotted contours), and SNe distance data (magenta dotted con- 
tours). The two contour levels in each case enclose regions con- 
taining 68.27% and 95.45% of the total likelihood. 



function as 

r{z) = 



Hiz'] 



C (zi,,, 



dz' 



H{z,) 



+ 



H(z„) 



(3) 



Expressing this relation as a linear interpolation rather than 
a stepwise function makes little difference to the results. The 
availability of the WiggleZ and BOSS Alcock-Paczynski dis- 
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Figure 7. Panels displaying, from left to right, measurements of the Hubble parameter H(z), the cosmic expansion rate a = H{z) /(I + 2), 
and the "Om" statistic [[H{z)/Ho]'^ - l] / [(1 + - l] fit as a stepwise function in 9 redshift bins of width Az = 0.1 using a Monte 
Carlo Markov Chain. The fit is performed to the WiggleZ and BOSS joint BAO and AP datasets, other BAO measurements from 6dFGS 
and SDSS, and SNe distance data. The measurement in each bin is marginalized over H(z) in the other bins, the spatial curvature f2]j 
and a WMAP prior for the sound horizon at baryon drag. The solid lines are not fits to the data but represent a fiducial flat ACDM 
cosmological model with parameters Om = 0.27 and h = 0.71. We do not plot a value for Om(2) in the first redshift bin because the 
statistic is not well-defined in the limit 2 — >■ 0. 



tortion measurements, with their direct sensitivity to H{z), 
brings two significant benefits to this analysis: a more pre- 
cise determination of H{z) in stepwise bins, and a lower 
covariance between the measurements in different bins. 

The left-hand panel of Figure 7 illustrates the measure- 
ments of H{z) m N = Q stepwise redshift bins of width 
A« = 0.1, where the likelihood is computed using the Wig- 
gleZ and BOSS joint BAO and AP datasets, the other BAO 
measurements from 6dFGS and SDSS, and the SNe Union 
2 dataset. We note that we do not use the WMAP distance 
priors in this fit because of the uncertainty in extrapolating 
the expansion rate beyond our bins, in the redshift range 
< z < Zt, m order to deduce the value of Da{z,) 
which is required to evaluate the quantities Ia and TZ. How- 
ever, we do marginalize over a WMAP-motivated Gaussian 
prior in Q.mh? with mean 0.1345 and width 0.0055. This 
prior is used when fitting to the BAO dataset, both when 
using the acoustic parameter A{z) oc vHrn/i^ and when cali- 
brating the sound horizon scale rs{zd)- In the latter case we 
assume a fiducial baryon density Q,\,h? = 0.0226. We also 
marginalize over a wide uniform prior in spatial curvature 
-1< fik < 1. 

We obtain measurements of H[z) with precision better 
than 7% in most = 0.1 redshift bins in the range z < 0.8 
(we note that the improved accuracy in the 0.7 < z < 0.8 bin 
compared to the adjacent bins is due to the presence of the 
WiggleZ Alcock-Paczynski data point at z = 0.73). As dis- 
played in Figure 7, our measurements are consistent with a 
fiat ACDM cosmological model with parameters Sim = 0.27 
and h = 0.71 (the value of the chi-squared statistic cal- 
culated using the full covariance matrix is 7.52 for 9 de- 
grees of freedom). Figure 8 displays the covariance between 
the measurements of H{z) in each bin, deduced from the 
MCMC chain. The correlation coefficients between different 
bins vary depending on whether AP data is available, but 
are generally low or moderate, r < 0.5. 

Table 3 lists the marginalized measurements of H{z) in 
each bin. We also convert these measurements to values of 
the cosmic expansion rate a = da/dt in physical units, where 
a = 1/(1 + z) is the cosmic scale factor, and values of the 



"Om" statistic (Sahni, Shafieloo . 
is defined by 



Starobinsky 2008) which 



Om(z) 



[H(z)/Hof - 1 
(1 + ^)3-1 ■ 



(4) 



In a spatially-fiat ACDM model this statistic is constant at 
different redshifts and equal to today's value of the matter 
density parameter flm. In universes with different curvature, 
or containing dark energy with different properties to a cos- 
mological constant, Om(2:) would evolve with redshift. These 
determinations of d and Om(2:) are plotted as the central and 
right-hand panels in Figure 7, respectively. We see a signif- 
icant decrease in the value of a between redshifts 2 = and 
= 0.7, corresponding to accelerating cosmic expansion. For 
example, the low-redshift expansion rate a = 0.069 ± 0.002 
Gyr"^ has dropped to d = 0.060 ± 0.002 Gyr~^ aX z = 0.55 
and d = 0.058 ± 0.004 at z — 0.75. The measurements of 
Om{z) are consistent with a constant ~ 0.25, as expected in 
a spatially-flat ACDM model. 

In Figure 9 we compare the results of fitting H{z) in 
5 bins of width Az = 0.2 to different subsets of the total 
dataset. We note that, when combined with a CMB prior, 
the baryon oscillation scale is calibrated in units of Mpc 
and permits direct measurement of H{z) in units of km 

Mpc~^. However, given that the normalization of the 
supernova Hubble diagram is treated as an unknown pa- 
rameter, the SNe data yield the relative luminosity distance 
Dl{z)Ho/c and require an extra prior in Ho in order to de- 
termine the function H{z). We take this prior as the Riess 
et al. (2011) 3% determination of the Hubble constant using 
new observations of Cepheid variables combined with Type 
la supernovae and the megamaser host galaxy NGC 4258, 
which yields a Gaussian prior in Ho of mean 73.8 km s~^ 
Mpc~^ and standard deviation 2.4 km s~^ Mpc~^. Figure 9 
illustrates the factor 2-3 gain in precision at higher redshifts 
achieved when adding the WiggleZ and BOSS AP measure- 
ments, with their direct dependence on H{z), to the existing 
BAO and SNe dataset. It is also notable that the differ- 
ent probes produce consistent determinations of the expan- 
sion history within the statistical errors, which are consis- 
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Table 3. Measurements of the Hubble parameter H{z), cos- 
mic expansion rate a = H(z)/{1 + z) and "Om" statistic 
[[H{z)/Ho]'^ - l] / [(1 + 2)^ - 1] fit as a stepwise function in 9 
redshift bins of width Az = 0.1 using a Monte Carlo Markov 
Chain. The fit is performed to the WiggleZ and BOSS joint BAO 
and AP datasets, other BAO measurements from 6dFGS and 
SDSS-DR7, and SNe distance data. The measurement in each 
bin is marginalized over H(z) in the other bins, the spatial cur- 
vature Qj; and a WMAP prior for the sound horizon at baryon 
drag. We do not quote a value for Om{z) in the first redshift bin 
because the statistic is not well-defined in the limit 2 — ^ 0. 



2 


H{z) 


a{z) 


Om{2) 




[km s~i Mpc^i] 


[Gyr-i] 




0.05 


70.2 ± 1.7 


0.0685 ±0.0017 




0.15 


77.4 ±2.7 


0.0687 ± 0.0024 


0.34 ±0.15 


0.25 


76.3 ±4.7 


0.0624 ± 0.0038 


0.15 ±0.14 


0.35 


81.9 ±4.5 


0.0620 ± 0.0034 


0.22 ±0.10 


0.45 


88.4 ±5.5 


0.0624 ± 0.0039 


0.26 ± 0.09 


0.55 


90.9 ±3.4 


0.0600 ± 0.0022 


0.23 ±0.04 


0.65 


103.7 ± 25.2 


0.0643 ±0.0156 


0.32 ±0.29 


0.75 


98.2 ± 6.2 


0.0575 ±0.0036 


0.21 ± 0.05 


0.85 


136.8 ±28.5 


0.0758 ±0.0158 


0.50 ±0.28 




H(z,) H(Z2) H(Z3) H(z,) H(Z5) H(z^) U(z^) H(2g) H(Zg) 



Figure 8. The matrix of correlation coefficients for the measure- 
ments of the Hubble parameter H{z) in 9 redshift bins plotted in 
Figure 7, obtained using a Monte Carlo Markov Chain. 



tent with a flat ACDM cosmological model with parameters 
r2m = 0.27 and h = 0.71. 

Another benefit of using the AP dataset is a reduced co- 
variance between the H{z) measurements in different bins, 
which are heavily correlated when only total distance infor- 
mation is available. This is illustrated by Figure 10, which 
displays the covariance between the H{z) measurements in 
5 bins for the case of SNe data + Ho prior, illustrating the 
significantly higher correlation coefficients r > 0.5 in com- 
parison with Figure 8. 

In order to illustrate further the consistent results ob- 
tained when varying the input datasets we fitted each de- 
termination of H{z) plotted in Figure 9 for the normaliza- 
tion factor h, assuming Q^n = 0.27, using the correspond- 
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Figure 9. Measurements of the Hubble parameter H(z) fit as 
a stepwise function in 5 redshift bins of width Az = 0.2 using 
a Monte Carlo Markov Chain. Results are compared for datasets 
SNe + Ho prior, SNe + AP + Ho prior, BAO + Vs (za) prior, BAO 
+ AP + rg{zii) prior. The measurement in each bin is marginal- 
ized over H{z) in the other bins and the spatial curvature fiij. 
The solid line is not a fit to the data but represents a fiducial 
fiat ACDM cosmological model with parameters = 0.27 and 
h = 0.71. 




H(z,) HCzj) H(Z3) H(2^) HCzj) 



Figure 10. The matrix of correlation coefficients for the mea- 
surements of the Hubble parameter H(z) in 5 redshift bins using 
just the SNe dataset and Ho prior, obtained using a Monte Carlo 
Markov Chain. 

ing covariance matrix obtained from the Markov chain. 
The resulting measurements for the cases (SNe, SNe+AP, 
BAO, BAO+AP) were h = (0.72 ± 0.03, 0.72 ± 0.03, 0.69 ± 
0.02, 0.68 ±0.02), respectively. [We note that the aim here is 
not to measure h, but to demonstrate that analyzing these 
subsets of the data produces consistent results.] 

4.4 Kinematical model fits 

Finally, we fitted our dataset with a "kinematical" cosmo- 
logical model (Rapetti et al. 2007) which is parameterized 
in terms of the dimensionless second and third derivatives 
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of the scale factor a(t) with respect to time, the decel- 
eration parameter q{t) — —H~^(a/a) and jerk parameter 
j{t) = H~^{'a /a). In particular, we adopt a parameteriza- 
tion where models are expressed in terms of the present-day 
value of the deceleration parameter qo and a constant jerk 
j, noting that ACDM models correspond to the special case 
j = 1. Following Rapetti et al. (2007), for a given {qo,j) we 
determine the function 



V{a) 



P - 



2p 



p + u 

~2r 



(5) 



where p = + 8j and u — 2{qo + \). Given the function 
V{a) we can determine the expansion rate as 

(6) 



Hiz) 



Ho 



2V{a) 



We note that there is a region in the (go, j) plane defined by 

j<qo + 2q^ qo < -1/4 

j<-l/8 qo >-l/4 (7) 

for which the condition V{a) > is not satisfied for all a 
and hence there is no Big Bang in the past. We exclude this 
region from our fits. 

Figure 11 illustrates the joint likelihood of kinematical 
model fits to the BAO, AP and SNe datasets. As in previous 
Sections, we marginalize over a WMAP-inspired Gaussian 
prior in Q^nh^ in order to calibrate the baryon oscillation 
standard ruler. ACDM models correspond to the line j — 1, 
and specific values of flm pick out a point with qo — l^m — 1. 
These models are consistent with the data. We note that 
SNe provide the best constraints on the kinematical model 
parameters, and that including the WMAP distance priors 
in the fitted dataset produces a very accurate joint con- 
straint on {qo,j) from the precisely-known distance to the 
last-scattering surface, although this corresponds to a signif- 
icant extrapolation of the validity of the model from a > 0.5 
to a > 0.001. Fitting to the combination of SNe, BAO and 
AP data, not including the CMB, produces marginalized pa- 
rameter measurements qo — — 0.67±0.16 and j = 1.37±0.68. 



5 CONCLUSIONS 

We have used large-scale structure measurements from the 
WiggleZ Dark Energy Survey to perform joint fits for the 
baryon-oscillation distance scale quantified by the acoustic 
parameter A{z) cx [Da{z)^/H(z)Y^^, the Alcock-Paczynski 
distortion parameter F{z) oc Da{z)H{z) and the normalized 
growth rate / erg (z) in three overlapping redshift slices with 
effective redshifts z = (0.44,0.6,0.73). We use lognormal 
realizations to quantify the covariances between parameters 
and redshift slices, producing a 9 x 9 covariance matrix. 

By combining the joint measurements of A{z) and F{z) 
taking into account the covariance, we performed simultane- 
ous determinations of the angular-diameter distance Da{z) 
and Hubble parameter H{z) based only on the WiggleZ Sur- 
vey dataset and a WMAP prior in the matter density il^nh? 
to calibrate the baryon oscillation standard ruler. These re- 
sults are generally insensitive to the fiducial cosmological 
model including spatial curvature. We measure these pa- 
rameters with 7 — 9% accuracy in each redshift bin. 

We use a combined dataset consisting of these joint 
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Figure 11. The joint probability of kinematical model param- 
eters go and j fitted to different combinations of datasets and 
assuming Hj; = 0. Results are shown for SNe data alone, the 
BAO -I- AP dataset and BAO -I- WMAP. The two contour levels 
in each case enclose regions containing 68.27% and 95.45% of the 
total likelihood. 



WiggleZ geometric measurements, other BAO data and SNe 
luminosity distances to perform a Monte Carlo Markov 
Chain determination of the expansion history H{z) as a 
stepwise function in 9 redshift bins of width Az — 0.1, 
also marginalizing over spatial curvature Slk. The results 
are consistent with a flat ACDM cosmological model with 
parameters Qm = 0.27 and h = 0.71. The addition of 
the AP data reduces both the errors in these measure- 
ments [through its direct sensitivity to H{z)] and the co- 
variance between different redshift bins. When we convert 
our results to a measurement of the cosmic expansion rate 
d = H{z)/{1 + z), we see a significant decrease in the value 
of a between redshifts z = and — 0.7, corresponding to 
accelerating cosmic expansion. Measurements of the statis- 
tic Om{z) = [[H{z)/Ho]'^ -l] / [(1 -I- zf - l] are constant 
with redshift, consistent with a spatially-flat ACDM model 
with matter density parameter fim « 0.25. 

We compare our measurements to cosmological mod- 
els including different dark energy equations-of-state w and 
kinematical models expressed in terms of the derivatives of 
the cosmic scale factor, the deceleration and jerk parame- 
ters. We find all data to be consistent with a cosmological 
constant model. 
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